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ABSTRACT 

The existence of bimodal disks is investigated. Following a simple argument based on energetic consid- 
erations we show that stationary, bimodal accretion disk models in which a Shakura-Sunyaev disk (SSD) 
at large radii matches an advection dominated accretion flow (ADAF) at smaller radii are never possible 
using the standard slim disk approach, unless some extra energy flux is present. The same argument, 
however, predicts the possibility of a transition from an outer Shapiro-Lightman-Eardley (SLE) disk to 
an ADAF, and from a SLE disk to a SSD. Both types of solutions have been found. We show that for any 
given accretion rate a whole family of bimodal SLE- ADAF and SLE-SSD solutions can be constructed, 
each model being characterized by the location of the transition radius. 

Subject headings: accretion, accretion disks — black hole physics — hydrodynamics 



1. INTRODUCTION 

Observations of black hole X-ray binaries (BHXBs, 
e.g. V404 Cyg, Nova Muscae, A0620-00, GRO J1655- 
40, Cyg X-l) suggest that a bimodal accretion disk 



may be present in these sour ces ( Thorne & Price 1975 



Shap iro, L ightman & Eardley 1976 ; Nar ayan, McClintock 
& Yi |l996fc Esin, McClintock & Narayan p97| ) . The inner 
part is a hot, optically thin, quasi-spherical accretion flow, 
producing a non-thermal spectrum via synchrotron radia- 
tion and Comptonization, with a cut-off around 100 kcV. 
The outer part is a ge ometr ically thin, optically thic k disk 
(Shakura & Sunyaev |l973| ; Frank, King & Raine |1992[) 
which emits a multi-color blackbody spectrum peaked 
around a few keV. By making plausible assumptions on 
how the transition radius Rt r depends on the accretion 
rate M, one obtains a one-parameter family of spectra, 
once the black hole mass M and the viscosity parame- 
ter a are specified. The different spectral states of the 
BHXBs can then be explained by varying the accretion 
rate M = rhMEddi where MEdd is the Eddington accre- 
tion rate. 

Historically the hot, inner region was first thought to 
be a cooling-dominated, optically thin a ccreti on disk, the 
SLE disk (Shapiro, Lightman & Eardley 1976 , henceforth 
SLE). However, such disks are subject to a violent ther- 
mal instability, which makes them useless for constructing 
stationary bimodal models. An alternative to the SLE 
disk could be the recently discovered hot, adve ction domi 



disk (SSD) connects smoothly to the inner ADAF, seems 
very promising in explaining the whole range of spectral 
states of BHXBs, from the quiesc ent to the high state 
(Esin, McClintock & Narayan |l997| ) . 

Advection dominated accretion flows are dynamically 
very stable. The robustness of the ADAF solutions led 



Narayan & Yi (|1995|) to suggest that perhaps nature will 
always select this option whenever it becomes available. 
Yet, the precise mechanism through which the SSD ma- 
terial is converted into an ADAF remains a matter of 
debate. Observations of Cyg X-l in the low state sug- 
gested that the transition from a cold to a hot disk 
might result from the secular instability of the radiation- 
pressure dominated inner part of the SSD (SLE; Ichimaru 
1977). The even more violent thermal instability of the 



radiation -press ure do minate d region (Pringle, Rees & Pa- 
cholczyk 1973 ; Piran 1978| ) might also be held respon- 
sible for such a transition. The drawback is that the 
transition radius for these models is usually quite close 



to the black hole, Rt 



45a^ 21 m 16 / 21 ml /21 R 



where 



M/Mq and R g = 2GM /c 2 is th e Schwarzschild 



radius (SLE; Frank, King & Raine 1992), whereas obser 



nated accretion flo ws (AD AFs: Narayan & Yi [1994 [1995 
Abramowicz et al. 1995 ), which are stable against both 
thermal and viscous perturbations. ADAF models have 
been quite successful in reproducing the spectral proper- 
ties of BHXBs in quiescence and of the G alacti c center 
source Sgr A* (Narayan, Yi & Mahadevan 1995 ). A bi- 
modal disk structure, in which an outer Shakura-Sunyaev 



vati ons of many BHXBs ( e.g. N arayan, Barret & McClin- 
tock 1997] Hameury et al. 1997) seem to imply a transition 
radius R tr ~ 10 4 R g . Moreover, no self-consistent global 
solution for such a viscously/thermally driven transition 
has been found so far, and, more probably, the thermal 
instability results in a genuinely time-dependent be havior 
(e.g. a limit cycle, see Mey er 
Mincshigc 



Wheeler 1989 



Meyer-Hofmeister 1981 
Honma 



Matsumoto & Kato 
|1993 ; Szuszkiewicz & Miller 1997a] , |j), rather than in a 
stationary transition to a hot flow. More promising in 
this respect are the evaporation models in which matter 
evaporates from the SSD forming an advec tion dominated 



corona (Meyer & Meyer-Hofmeister 1994; Narayan & Yi 
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1995) 



The question remains whether we understand why it 
seems not to be possible to convert an SSD directly into 
an ADAF at a certain radius. This is the question we ad- 
dress in this paper. We use a simple argument based on 
energetics to predict which kind of bimodal disk structures 
should be possible. We conclude that only bimodal accre- 
tion flows which are formed by a SLE disk outside and by 
an ADAF or SSD inside are allowed. We confirm our qual- 
itative prediction by actually constructing global models 
for both bimodal SLE- ADAF and SLE-SSD solutions. 

2. THE EQUATIONS 

The system of equations governing the disk structure 
used in the present investi gation are basically those of slim 
disks (Abramowicz et al. 1988) and our appro ach fo llow 
closely those used by Narayan, Kato & Honma ( [1997 ) and 
Chen, Abramowicz & Lasota (1997) in investigating the 
properties of global ADAF solutions. More specifically, 
we assume that the disk is stationary, axisymmetric, that 
the gravitational field of the central object is described 
by the pseudo-Newt onian potential $ = —GM /(R— R g ) 
(Paczyhski & Wiita 1980] ) and that self-gravity is com- 
pletely negligible. Contrary to Narayan, Kato & Honma, 
we include radiative cooling in our equations, using the 
approximate expre ssions for bremsstrahlung losses given 
by Narayan & Yi (199!:), which are valid at any optical 
depth. Although we expect optically thin bremsstrahlung 
cooling and gas pressure support to be adequate for con- 
sistently modeling the SLE- ADAF, the Shakura-Sunyaev 
disk in SLE-SSD models may be optically thick and radia- 
tion pressure dominated, so we include radiation pressure 
in our calculations. 

By denoting the total pressure and the mass density 
at the equatorial plane with p and p, the corresponding 
vertically-integrated quantities are given by 



P = 2Hp , E = 2Hp 



(1) 



where H is the half-thickness of the disk. The pressure p 
is the sum of gas and radiation pressure 



P = Pg+Pr, 

K 

Pg = -pT 



(2) 
(3) 



where T is the average temperature of electrons and ions 



T = fi[— + — 

,Pi Pe 



(4) 



and for simplicity we take /ij = p, e = 1, p, = 0.5. We 
allow Ti and T e to deviate, but instead of using a so- 
phisticated model for electron-ion coupling we just as- 
sume T e = min(Ti, 6 x 10 9 K). Following Abramowicz 
et al. ( 1996 ), we adopt for the radiation pressure the ex- 
pression 



Pr 



Q- 

4c 



2 

V3 



(5) 



where Q- is the cooling rate and r is the total optical 
depth (see below). Although in principle also magnetic 
pressure should be accounted for, this effect, like all other 



effects of the frozen-in i?-field, will be ignored in this pa- 
per. The isothermal sound speed c s is given by 



r 2 _ P P 

s S p 



(0) 



and, following the standard approach, the disk thickness 
and the kinematic viscosity are expressed as 



H = 



K 



v = -ac s H 



(7) 



(8) 



where VLk is the Keplerian rotational frequency, calculated 
using the pseudo-Newtonian potential. Introducing the 
radial velocity vr and the angular velocity f2, the continu- 
ity, radial momentum, azimuthal momentum and energy 
equations take the form 



vr 



ZRv R 



Al = -2-kTjRvr . 



(n 2 -n 2 K )R 



dR 

d{R 2 n) _ d 

~ dR 
dp- 1 



1 dp 
~pdR 



0. 



dR 



Y<vR 



.dtt\ 



p- 



dR 



dR) ' 
Q+-Q-. 



(9) 
(10) 

(11) 
(12) 



We adopt for the specific internal energy (including radia- 
tion) an equation of state of the form (7 — l)e — p/ p = c 2 
with constant adiabatic index. In the energy equation 
(eq. [p2| ) , the viscous heating Q + has the usual expression 



(13) 



while the cooling rate Q_ is calculated using a bridging 
formula valid both in th e opti cally thin and optically thick 
regimes (Narayan & Yi 1995| ), 



Q- = 8aT* 



St 



V3- 



8(rTi 



(14) 



2 ' ' ~ ' Qt r 
Here r is the total (scattering plus absorption) depth 

T = T es +T abs (15) 

and 

Tes = ^esS , (16) 

Q b l r 



Tabs — 



8<tT 4 



(17) 



The bremsstra hlung cooling function is expressed as 
(Narayan & Yi |199E 



(18) 



Q b I = 2Hq b l = 8.5 x lO^ff^/e^ergcm-V 1 
where 

f ei (6 e ) = 4 f^V (1 + 1-781 e 34 ) (19) 



and 9 e 
aturc. 



kT e /m e c 2 is the dimensionless electron temper- 
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3. BIMODAL DISKS: ENERGETICS 

The structure of accretion disks is determined by the set 
of differential equations given in section ||. In the case of 
geometrically thin disks these equations become algebraic, 
so that the disk properti es ma y be evaluated locally (see 
e.g. Frank, King & Raine p.992| ). Although this is not true 
any more for geometrically thick disks, a good estimate 
of their structure may be nevertheless obt ained from local 
considerations (see e.g. Abramowicz et al. 1995| ). 

In the analysis of the global behaviour of the models, 
a key role is played by the energy equation, which simply 
states that viscous heating, radiative cooling and radial 
energy advection must balance at any radius 



= Q+ — Q- — Qadv — . 



(20) 



The non-locality of equation ( J20j ) is contained in both 
Qadv and Q+, which are functions of the gradients. For 
a standard thin disk, however, Q a dv — and Q+ oc rh, 
while for the other classes of disks they can be reasonably 
well approximated in terms of local variables. With this 
proviso, the energy equation, at any given radius, may be 
written as a function of rh and S: Q = Q(S,rh). The en- 
ergetically permitted configurations are then given by the 
intersection of the surface Q(E,rh) with the Q = plane 
as shown i n figu re [I] (see e.g. Abramowicz et al. 191 
Chen et al. 1995). The curve on the left corresponds to 
ADAFs (upper branch) and SLE disks (lower branch) , the 
S-shaped curve to 'slim disks' (upper branch), radiation 
pressure supported SSD (middle branch) and gas pressure 
supported SSD (lower branch). 

Two of the five branches of accretion disks are well 
known to be thermally unstable, as can be easily inferred 
from figure HI. The unstable branches are those for which 
Q + > Q- above and Q + < Q_ below the curve. For a 
solution on the radiation pressure supported SSD or SLE 
branch, a sudden local heating (which amounts to a ver- 
tical upward displacement in the log £ — log rh diagram) 
causes the disk to enter a region where heating exceeds 
cooling (Q + > Q-). This results in more heating and a 
runaway thermal instability sets is. A local cooling of the 
unstable disk similarly produces a runaway cooling. Al- 
though the plot of figure 1 refers to the low viscosity case, 
the same conclusions hold for large values of a ~ 0.2 — 1, 
when the topology of the solutions in the log £ — log rh 
plane changes (a very clear example can be found e.g. in 
Abramowicz & Lasota 1995). 

In the rest of this section we show how the log £ — log rh 
diagram can be used also to identify the possible classes 
of bimodal, stationary flows, i.e. two types of disks con- 
nected by a boundary layer. Although the following anal- 
ysis is somewhat qualitative, it leads to rather interesting 
conclusions which are confirmed by the numerical calcula- 
tions discussed in sections [| and |[ Let us consider again 
the energy equation. Using the exact (differential) expres- 
sion for Qadv, but retaining a 'local' approximation for 
Q+, equation (p0|) becomes 



ladv = SWflT— = Q + - Q_ 

an 



(21) 



where s is the specific entropy. In a stationary flow 
vnds/dR is just the Lagrangian rate of change of the spe- 



cific entropy of a gas element (in the case under examina- 
tion, a 'vertically-integrated' thin ring). Suppose now that 
the gas parcel is located at a given radius and that it has 
a specified value of s = Sq. Assume, just for concreteness, 
that so is such that the configuration lies a bit to the left of 
the stable SSD branch of figure |l], so that it is Q + < Q_. 
As time elapses, viscous stresses cause the gas parcel to 
move inwards and, since rh is constant, this amounts to 
a horizontal shift in the log £ — log rh diagram. Because 
cooling exceeds heating, the specific entropy must decrease 
in time and this implies that the configuration moves to 
the right, i.e. towards the curve. If the element is initially 
to the right of the branch, heating exceeds cooling, s in- 
creases in time and the gas parcel shifts again towards the 
SSD branch. One can then conclude that, even if the flow 
is not in the SSD phase at a certain radius, it will try to 
become a SSD for smaller radii. The same reasoning holds 
both for the stable and for the unstable SSD branches, 
since, in both cases, one has Q+ < Q_ to the left and 
Q+ > Q- to the right of the curve. A major consequence 
is that stationary, bimodal accretion flows which consists 
of an outer SSD and a different kind of flow at smaller 
radii seem not to be allowed. 

When applied to the SLE branch, however, this argu- 
ment yields just the opposite result. Now it is Q + < Q- 
to the right and Q+ > Q- to the left of the curve, so, 
unless it lies exactly on the SLE branch, the spiraling-in 
matter moves away from the curve. Depending on the ini- 
tial value of the entropy, it will end either on the ADAF 
or on the SSD branch. Bimodal, stationary SLE-ADAF 
and SLE-SSD models seem therefore possible, in the sense 
that they are at least energetically permitted. 

The same conclusions may be drawn looking at the in- 
tegral curves of equation (|2ll) , which are shown in figure 
|3[ By integrating equation (|2l|) inwards, starting close to 
the SLE branch, the solutions diverge exponentially, until 
either the ADAF or the SSD branch is reached. This is the 
stationary transition from a SLE to either an ADAF or a 
SSD. A similar transition cannot be found starting from 
a SSD/ ADAF and integrating inwards, because now the 
solutions exponentially relax back onto the SSD/ ADAF 
branch. 

In general one can say that disks on a branch with 
Q+ > Q- on its left and Q + < Q_ on its right can per- 
form a stationary transition to a different state towards 
smaller radii, while disks on a branch with Q- > Q+ on 
its left and Q_ < Q + on its right can not, as illustrated 
in figure ^. Much as in the case of the stability analy- 
sis (see above), the same considerations apply to high-a 
flows (as can be eas ily verified from figure 1 of Abramow- 
icz & Lasota 1995 ), so our result is independent of the 
viscosity parameter. Clearly, the inclusion of some extra 
(non -local) heating/cooling process in the energy equation 
chang es the picture outlined above. As show n by Honma 
( 199q ) and Meyer & Meyer-Hofmeister (1994), convection 
and electron conduction may indeed provide a way to ob- 
tain global SSD-ADAF models. 

The analysis presented in this section relies on the tacit 
assumption that none of the branches change as we follow 
the evolution of a gas element towards smaller radii. Of 
course this is not entirely true, since, at variance with what 
is done in the classical stability analysis, the radial coor- 
dinate of the gas parcel now varies because matter is spi- 
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raling in. Strictly speaking, one should draw many curves 
on the log E — log rh diagram, each referring to the differ- 
ent radii crossed by the gas element. We note, however, 
that the time over which the topology of the log E — log m 
diagram changes is the viscous time of radial motion. The 
gas element evolves on a thermal timescale which is much 
shorter than the viscous time scale, so the small defor- 
mation of the curves can be safely ignored to first order. 
Much more delicate is the assumption that Q+ can be eval- 
uated in terms of local variables. In the transition region 
the derivative of Q may deviate strongly from its usual 
value, casting some doubts on our conclusions. A strong 
test for the simple argument presented above is to actu- 
ally construct the bimodal solutions for the cases in which 
transitions are allowed, by solving the full set of differen- 
tial equations. Both the global SLE-ADAF and SLE-SSD 
solutions have been found numerically and are presented 
in sections | and |. 

4. BIMODAL SLE-ADAF FLOWS 

In this section we present our numerical models for bi- 
modal SLE-ADAF disks. We construct solutions which 
are in the ADAF phase for small radii and in the SLE 
phase for large radii. We shall show that the radius Rt r at 
which the transition occurs is a free parameter, with the 
only limitation that it must be Rt r < Rmax ■ The existence 
of a maximum transition radius follows from the fact that, 
for too large radii, radiative cooling becomes too strong to 
allow of a hot, optically thin solution. A rough estimate 
of Rmax may found by equating the viscous heating to the 
optically thin bremsstrahlung cooling (eqs. |l3| and JLst) 
and giving reasonable guesses of the values ofvgjR), Cs(R) 
and for an ADAF (Abramowicz et al. |1995|), 



R-n 



2 x lO^R g . 



(22) 



Sufficiently within this radius bimodal stationary solu- 
tions, which are SLE for R tr < R < R m ax and ADAF 
for R < R tr , exist. 

Equations (|l^, [ll], |l^) form a fourth order set of dif- 
ferential equations. We write the second order angular 
momentum equation ( |TT| ) as two first order equations and 
solve the resulting four coupled differential equations nu- 
merically using a Henyey-type relaxation scheme (Nobili 
& Turolla |L988|) . Since we do not use the first integral of 
equation (|11|), the specific angular momentum at the in- 
ner edge, £q, does not appear exp licitly as an eigenvalue 
(see e.g. Narayan, Kato & Honma 1997), but follows from 
the solution once suitable boundary conditions have been 
specified. Logarithmic variables are adopted throughout: 
lnr, where r = 2R/R g , ln(vu/c), ln(c s /c), ln(£l/f2^) and 
dhiQ/dhir. 

Global ADA F mod els are transsonic (see e.g. Na rayan , 
Kato & Honma 1997; Chen, Abramowicz & Lasota 1997) 
and we expect also the bimodal solutions we are looking 
for to pass through a critical point close enough to in- 
ner rim. The location of the sonic point is an eigenvalue 
of the problem and follows once all boundary conditions 
are specified. The Henyey solver we use is specifically de- 
signed to locate critical points automatically, looking for 
the points where the Jacobian of the system changes sign. 
When a critical point is met, the solver replaces one of the 



boundary condition with a regularity condition. Like in 
all relaxation-type schemes, a trial solution must be sup- 
plied to start the iteration. Usually its exact form is not 
very important, the only strong requirement being that 
the initial guess is transsonic. In the present case, we first 
generated a global ADAF model, using a modified self- 
similar solution as first guess, and then used it as the trial 
solution for the bimodal disk. 

Although different values of the disk parameters were 
tested, models discussed here were obtained for a — 0.1, 
7 = 1.5, rh — 10~ 3 and m* = 10. The inner and 
outer edges of the disk are located at Ri n — l.5R g and 
Rend = l0 6 R g , respectively. An equally spaced logarith- 
mic radial grid was used with N ~ 300 points. Grid- 
refinement was used to obtain better spatial resolution in 
the transition region. 

4.1. The Boundary Conditions 

In order to clarify the choice of the boundary condi- 
tions, let us start discussing how boundary conditions 
should be placed for a global ADAF solution in our scheme. 
Since the structure equations translate into four first or- 
der ODEs, four conditions are needed. As discussed above, 
the Henyey solver automatically puts one of these condi- 
tions (a regularity condition) at the critical point, which 
is automatically located. If a fine enough grid is used, 
any constraint which forces the derivatives to stay finite 
at the sonic radius will do. The remaining three conditions 
can be placed at either R il } or R ou t- As discussed e.g. by 
Narayan, Kato & Honma (1997) the inner boundary con- 
dition must require the specific angular momentum, i, to 
behave physically. At Ri n , we ask that d£/dR = 0, which 
is equivalent to dlnil/dlni? = —2. The last two condi- 
tions are given at the outer edge, specifying the values of 
c s and f2 at R ou t (see again Chen, Abram owicz & Lasota 



1997| and Narayan, Kato & Honma |1997| ). 

For the bimodal solutions we investigate in this paper, 
the outer disk is of the SLE type. Such a disk allows 
only one boundary condition at the outer edge, and the 
remaining condition must be specified elsewhere. This is 
a somewhat subtle point which will be discussed below. 

4.2. The Solutions 

In order to compute the SLE-ADAF solutions we start 
from a global ADAF model with very small rh. At the 
outer edge the two 'standard' boundary conditions 



Q{Rout) — uQ,K{Rout) 
Cs(Rout) = oSIk (Rout) Rout 



(23) 
(24) 



are prescribed, where < a < 1 and < ui < 1 are two 
arbitrary dimensionless constants. In the following, only 
the case uj = 1 is considered for the sake of simplicity. 
Once these global ADAF solutions have been computed, 
rh is gradually increased up to ~ 10 -3 , each time using 
the previous model as the input to the next run. For 
Rout — l0 6 R g this does still not influence the solutions 
significantly, provided that a remains close to 1. Starting 
from these solutions, models with smaller a are calculated. 
We find, however, that, owing to the increasing importance 
of bremsstrahalung cooling, no solution exists below a cer- 
tain value of a, which depends on rh and R ou t- As soon as 
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this happens, we move the outer boundary condition for 
c s (eq. p4|| ) to a smaller radius R tr < R ou t, 



C a (Rtr) = <jfl K (Rtr)Rt 



(25) 



using at the same time a value of a just above 'critical' 
[a ~ 0.3 typically) . One of the two outer boundary condi- 
tions (eq. |24|) is therefore transformed into an 'internal' 
condition (a constraint which is placed somewhere within 
the computational domain), which our solver can handle 
very easily. For any value of Rtr we find a global solution 
which turns out to be an ADAF for R < R tr and a SLE 
disk for R ^ R tr ■ We thus obtain a family of SLE-ADAF 
solutions parametrized by the value of the transition ra- 
dius Rtr- 

Our numerical solutions are shown in figure ^, where 
the disk height is plotted vs. radius for different values of 
R tr - Bimodal SLE-ADAF solutions have been fou nd pre - 
viously by Igumenshchev, Abramowicz & Novikov ( 1997| ), 
but, owing to a different choice of the boundary condi- 
tions, numerical problems prevented them from comput- 
ing the entire family of models. The reason is that all 
solutions are characterized by almost the same value of 
H = c s /£Ik at the outer edge and, at least those with 
large enough R tr , also at the inner edge, as can be seen 
from figure 0. This means that specifying c s (or H or any 
other variable) at the outer or the inner edge, instead of 
the internal condition used here, is formally correct but 
numerically unstable. Since the problem becomes nearly 
singular (many solutions with almost the same boundary 
condition), an impossibly accurate fine tuning would be 
required for many of the solutions we have found. By set- 
ting the condition at the transition radius this problem is 
solved. 

A close-up of the transition region is shown in figures ^[ 
|(], |?| and |[ Details in the behavior of the flow variables in 
the transition region are of importance inasmuch they can 
help in shedding light on the physical processes responsible 
for the transition itself. One of the peculiar features of the 
transition region is the small peak in the angular velocity 
which is seen in figure |^. For sufficiently small Rtr/ Rmax, 
as in the model considered here, super-Keplerian rotation 
may s et in. A similar behavior has been found by Honma 
( 1996| ) in modeling a SSD-ADAF transition with convec- 
tive fl ux, and by Igumenshchev, Abramowicz & Novikov 
(1997). Presumably the super-Keplerian rotation may be 
attributed to the pressure drop experienced by the gas as 
it ent ers the ADAF (Abramowicz, Igumenshchev & Lasota 
1997|). 



5. BIMODAL SLE-SSD FLOWS 

In this section we present numerical solutions for bi- 
modal SLE-SSD models. These are disks which are of 
SLE type for large radii and of Shakura-Sunyaev type for 
smaller radii. The choice of the parameters is the same as 
for the solutions discussed in the previous section, a = 0.1, 
7 = 1.5, m = 10~ 3 and = 10, although we also tried 
other values. As in the SLE-ADAF case the transition 
radius is arbitrary, as long as it is Rt r < Rmax- Now, 
however, we do not construct the full transsonic solutions, 
and take the inner radius sufficiently far away from the 
black hole that relativistic corrections are not essential for 
the flow behavior. 



The numerical approach is similar to that outlined in 
sectio n H , but now we use the angular momentum equa- 
tion (pT^ in its integrated form, 



vR 2 ^=v R (nR 2 -£ ). 



(26) 



This reduces the number of first order differential equa- 
tions from four to three. Since Ri n R g , the specific 
angular momentum at the inner edge can be assumed to 
vanish, £q — 0. We verified that any sufficiently small 
value is allowed, and found that the results are not very 
dependent on it. 

The SSD part of our bimodal SLE-SSD solutions is opti- 
cally thick and may be radiation pressure supported. Since 
the dynamical equations (|l^, O, 12) contains the total 
pressure, rather than the gas pressure, it is not straight- 
forward to compute the gas temperature T from a given 
value of p (or c s ). In order to get T, equations (||, [| |, |[ 
14 ) were solved in a separate iteration loop at each radius 
for a given value of the total pressure. 

5.1. The Boundary Conditions 



In contrast to what was discussed in section 4.1, the 



number of boundary conditions now is three and, since 
the solution is not transonic, none of them needs to be 
specified at a critical point this time. At the inner edge 
we want the flow to be a SSD, so we ask that the radial 
velocity there attains the value predicted by the standard 
disk model 



VR(Rin) = - 



](Rin) 



^K(Rin)Ri 



(27) 



As in the SLE-ADAF case, one boundary condition is 
given at the outer edge, where the flow is in the SLE phase, 
and, again, we take 



Q(Rout) — ^K(Rout) ■ 



(28) 



The remaining condition, the one on c s , needs to be speci- 
fied at the radius where we choose the transition to occur. 

5.2. The Solutions 

The way we actually obtained the bimodal SLE-SSD so- 
lution is similar to that described in section [|. We started 
out with a global SLE solution and then worked towards 
the bimodal solution by slowly modifying the inner bound- 
ary condition on the sound speed. By doing so we ran into 
the same problems as in the SLE-ADAF case, which could 
only be solved by moving the condition on c s to a point 
inside the computational domain, transforming it into an 
internal condition. The radius at which this condition is 
specified is again the transition radius. In the present case, 
however, the numerical integration is much more trouble- 
some because the equations become extremely stiff in the 
very thin transition region and in the SSD. This is not 
surprising since the cooling length in the SSD is extremely 
small, and in the transition region it is even smaller. To 
circumvent this problem we slightly modified the energy 
equation to guarantee an artificial lower limit to the cool- 
ing length. We believe that this does not invalidate our 
results, as long as the cooling length remains the small- 
est length scale of the problem. A particular solution is 
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shown in figure Contrary to the SLE-ADAF case, now 
it is very difficult to obtain solutions with a different value 
of Rtr, although possible. For this reason, and also be- 
cause of the little practical significance of these solutions, 
we just present here one model. 

6. DISCUSSION AND CONCLUSIONS 

The goal of this investigation has been to shed some 
light on the physics of disk transitions. Using a simple 
argument based on energetics, we have shown that global, 
stationary bimodal solutions are possible if the outer disk 
satisfies a simple condition. Our main result is that only 
an SLE disk can match an inner ADAF (or SSD), while a 
bimodal accretion configuration formed by an outer SSD 
and an inner ADAF is not permitted, at least within the 
standard slim disk physics. 

The two allowed bimodal structures, an outer SLE disk 
glued to an inner ADAF or SSD, have been computed nu- 
merically. We found that, in both cases, the transition 
radius turns out to be arbitrary, so a whole family of solu- 
tions exists. This may sound odd, since there is no physical 
reason to ask the sound speed to take a precise value at an 
intermediate radius, as we did in finding our numerical so- 
lutions. In principle the condition should be placed at the 
outer edge. In doing so, however, one should be extremely 
accurate in fixing c s (R ou t) to find the transition at the 
desired radius since most bimodal solutions correspond to 
very nearly identical values of c s (or any other variable) 
at the outer or inner edge. It remains nevertheless true 
that the structure of bimodal disks is determined once a 
complete set of conditions is prescribed at the edges and 



at the sonic radius, without the need of any extra piece of 
physics. 

The existence of the SLE-ADAF and SLE-SSD bimodal 
solutions agrees with the energy argument of section 3[ 
and strengthens its plausibility. The same argument indi- 
cates that a SSD-ADAF model is not allowed. There is a 
caveat here, however. The essence of the energy argument 
lies in recognizing that the non-locality of the advective 
cooling term in the energy equation can be removed by 
performing a transformation to the fluid frame, and doing 
the thermal analysis in the frame comoving with the fluid. 
If, on the other hand, thermal conduction is introduced, 
the energy equation becomes a second order diffusion type 
equation and there is no way of removing the non-locality 
by changing the observer's frame. In th is cas e our argu- 
ment definitely does not apply. Honma ( 1996| ) has shown 
that SSD-ADAF structures are possible when the thermal 
heat flux from the hot ADAF into the SSD is accounted 
for. In such accretion flows the location of the transition 
follows self-consistently from the model, contrary to what 
we have found in this paper. This is because the heat flux 
will evaporate per unit time an amount of matter from the 
SSD which needs to be in balance with the actual accretion 
rate. The transition radius will therefore automatically ad- 
just itself until this balance is reached. In concluding, we 
would like to state clearly that, because of the thermal in- 
stability of the SLE disk, the models presented here most 
probably do not exist in nature. 

We wish to thank I. Igumenshchev, M.A. Abramowicz 
and R. Narayan for useful comments and suggestions. 
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Log© 

Fig. 1. — Unified picture of accretion flows onto black holes. Solid lines mark the different branches of solutions. A branch is stable when 
arrows point towards the solid line. 




Fig. 2. — Same as in figure [j] for the transition analysis. Stationary, global bimodal models are possible for those solutions which lie on 
branches with outward pointing arrows. The only option is the SLE branch. The radiation pressure dominated SSD branch is unstable, but 
cannot change into a hot ADAF in a stationary model. 




Log(R) 

Fig. 3. — A schematic plot of the integral curves (dashed lines) of equation (|2l|). The solid lines are the ADAF, SLE and SSD branches. 
Each curve is found by starting at some radius with a given specific entropy ana integrating inwards. It is clearly seen that integral curves 
tend to move away from the SLE branch but converge to either the ADAF or the SSD branch. 
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Fig. 4. — The disk geometry for the family of SLE-ADAF models with m = 10 3 , m* 
characterized by a different location of the transition radius. 



10, a = 0.1 and 7 = 1.5. Each model is 




Fig. 5. — Close-up of the transition region for one of the SLE-ADAF transition solutions of figure |4|. The star marks the position of the 
transition radius. 




Fig. 6. — Same as in figure H for the specific angular momentum (solid line) and for the Keplerian angular momentum (dotted line). 
Super— Keplerian motion is clearly visible. 
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Fig. 7. — Same as in figure ]B for the density (in cgs units). 
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Fig. 8. — Same as in figure ^| for the pressure (in cgs units). The pressure drop as the SLE material enters the ADAF is clearly seen. The 
small wiggle at the critical point is a numerical artifact and is of no significance. 
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Fig. 9. — The isothermal sound speed (in units of c) for a bimodal SLE-SSD configuration. The star marks the position of the internal 
boundary condition Rtr- The solution is in the SLE phase for R > Rt r and in the SSD phase for R < Rtr- 



